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intensity  along  a raypath  in  terms  of  the  horizontal  dis- 


tance traveled  by  the  wave  front  and  the  angle  which  the 


ray  make a vith  the  horizontal.  Both  formulas  are  amenable 
to  digital  computation  in  parallel  with  a regular  digital 

iv.'  ath  calculation — 


For  over  a decade,  the  standard  cleans  for  computing  intensity  along  a 
sound,  raypath  in  waters  with  variable  sound  velocity  (versus  depth)  has 
been  to  approximate  the  sound  velocity  function  by  a continuous  broken-line 
function,  and,  then,  compute  the  vertical  distance  separating  two  closely 
adjacent  ray3.  'fhia  vertical  distance  is  compared  with  the  initial  angle 
between  one  rays  in  order  to  determine  an  approximate  value  of  the 
divergence . The  difficulties  with  this  approach  concern  the  sound  velocity 
appx’oximation  raid  the  closeness  with  which  the  two  neighboring  rays  oust 
be  taken  to  yStld  reasonable  intensities.  Obviously,  the  closer  a pair  of 
rays  aro  taken  in  initial  angle,  the  more  accurate  the  calculation  rsust  be 
for  each  ray  to  preserve  a significant  value  of  vertical  distance  between 
theu. 
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There  are  two  methods  by  which  the  standard  calculations  of  the 
pair jo?  rays  can  be  effected;  sequentially  or  simultaneously.  In  the 
former,  one  computes  the  first  ray  and  stores  the  2 (vertical)  oi'dinates 
for  a certain  predetermined  x (horizontal)  values;  then  the  same  is  done 
for  the  second  ray;  and  finally  the  stored  2 -ordinates  are  differenced 
to  fora  Az's.  In  the  latter,  one  has  two  computers  of  matched  accuracy 
for  the  computation  of  the  raypathn;  the  computation  yielding  the 
respective  z ordinates  for  the  same  running  variable  x;  then  a different 
ing  calculation,  either  continuous  or  discrete,  between  the  z’s  to  form 
Az  as  a function  of  x follows.  Both  means  require  high  accuracy  of  ray- 
path  calculation  (higher  than  if  one  Just  wished  to  compute  the  ray  it- 
self); both  ora  amenable  to  digital  calculation;  only  the  latter  seems 
suitable  for  analog  implementation.  The  former  requires  a very  extensive 
and  fast  memory  and  a large  computer  to  achieve  any  considerable  speed 
at  mailing  a complete  sound  raypath  and  intensity  plot.  The  latter  requires 
matched  accuracies  that  ere  difficult  and  expensive  to  obtain  and  maintain 
from  an  analog  ’unit. 

The  goal  behind  the  approach  initiated  here  is  to  find  an  auxiliary 
xeans  of  computing  the  sound  intensity  change  along  a raypath . The  fora 
of  the  calculation  would  be  to  calculate  the  coordinate’s  of  a single  my 
vfhilc  so  doing,  put  these  into  a second  equation  which  would  compute 
the  sound  intensity  (or  seme  function  of  it)  directly.  It  is  hoped  that 
in  this  manner , seme  of  the  high  accuracy  and  matched  unit  requirement a 
can  be  eliminated  and  reasonable  analog  means  arrived  at. 


Sound  Raypath  Intensity  Calculation 

Because  the  lateral  extent  of  the  oceans  is  many  times  greater  than 
the  depth  dimension,  the  fall  off  in  intensity  along  a ray  ia  measured 
in  terms  of  horizontal  aistanee  x.  Acoustic  intensity  is  power  per  unit 
area  and  its  loss  is  measured  from  some  standard  sphere  of  unit  radius. 

The  intensity  a,3  a fund  ic  _ of  x along  the  raynath  is  given  approximately 

by 


P coc  G (Aj/)(a9  ) 

I ?=  — — — — " 1 

x(A0)(Az)  cos  9 

where  the  ©,  x,  and  z quantities  are  as  shown  in  Figure  1 and  0 is  angle 
in  the  horizontal  plana.  As  can  be  seen,  the  usual  assumptions  apply: 
the  lirection  of  the  sound  raypath  does  not  vary  in  the  horizontal  plane, 
thuc  sound  velocity  is  n 2 dependent  function  only;  z is  ohe  depth 
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dimension  and  is  measured  with  increasing  size  downward;  0 is  the  angle 
the  raypath  makes  with  the  horizontal  axis  measured  positively  from  x 
to  z;  and  0 is  the  initial  angle  at  which  the  rey  leaves  the  trans- 
ducer. 


Depth  Plane  shoving  Gq,  0,  x,  z along  a Reypath 
Figure  1 

The  intensity  e Ion?  the  raypath  can  he  defined  c.n  follows: 


lim  P cos 
I 5 A0-*O 


AO  -X)  x cos 
o 


0 /AO  \ 

— ( — ) 
0 V Az  / 


P cos  0 


(?) 


(1) 


X COD  0 

where  it  oust  be  roceabered  to  tefce  Az  or  dz  in  the  positive  direction. 
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Eha  problem  of  evaluating  I retlucee  to  evaluating  d3  /d-.j  or 
inversely  dz/dO^. 

dz/iQ^  for  tae  Eg  Gradient  Case 

Oils  is  simply  done  f.  ui ref ersnce  to  the  ekatch . 


Az  =5  sec  0 (x  sec  S ) AO 
o c o 


2 

x sac 


0 

o 


dz 

Wo 


X 


3£C 


2 


0 

o 


Substituting  into  equation  (i)  one  gets 


P cos  0 


2 ~ 2 2 
x cos  0 xsec  © x sec  0 

00  o 


(2) 
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Equation  (2)  reduces  to  the  inverse  scuai'i  law  along  the  path  length, 
R,  when  one  notes  that 

R a x sec  0 . 

o 

dz/dOQ  for  the  Linear  Gradient  Case 


It  is  well  known  that  for  the  linear  gradient  case . the  raypaths 
become  arcs  of  circles . Let  the  velocity  equation  be  given  by 

7(z)  = VQ  + V (z  - zo)  (3) 

where  V is  the  3ound  velocity  at  depth  z^  and  V'  represents  the  gradient  — 
it  is  a°negative  number  for  decreasing  sound  velocity  with  increasing  depth. 

Snell's  law  rJLong  a particular  sound  ray  can  be  written 


But 


dz 

5x 


= tan  0 


and,  therefore.. 


dx  -V, 


cos  0. 


d©  V cos  0, 
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Differentiate  and  obtain 


or 


dz 


dz 

dO 


sin  0 


dO  cos  0 


-V 


sin  0. 


V cos  0 
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Let  us  define  Aq  by 


A_  S 


V'  cos  0 


(5) 


and  thus  we  get  as  the  diffential  equation  for  the  sound  raypath 
% = Ao  sin  0 


g = Ac  cos  0 


(6) 


which  are  evidently  circles  of  radius  A . Integrating  and  applying  the 
conditions  x = 0/  2 at  0 = 0q  givea°the  parametric  equations  of  the 
raypath  circle; 


z - x = -A  (cos  0 - cos  0 ) 
00  o 


1 

[ 


x =i  A (sin  0 - sin  0 ) 


(7) 


Elimination  of  the  parameter  0 gives  the  equation  of  the  circle  in  the 
closed  form 


0 * (z  - z )2  + x2  - 2A  cos  0,  (z  - z ) + 2A  sin  0 x.  (8) 
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Sound  Raypath  Circle 


Sound  Velocity  Curve 


It  is  evident  "ran  equations  (3)  and  vO  and  Figure  2 that  a physical 
raypath  will  not  travel  around  a complete  circle  because  to  do  so  requires 
a negative  velocity.  In  fact,  since  in  the  oce-uis  the  sound  velocity 
rarely  varies  by  more  than  8 o/o  total,  only  small  arcs  of  these  circles 
will  be  encountered  in  practice. 

The  computation  of  oiz/dO  is  done  as  a limit  of  the  ratio  Az/AOr 

Suppose  what  we  take  two  nearby  rays  eraitted  from  the  same  transducer,  a 

0 ray  and  a 0.  ray.  The  equation  for  the  circular  raypaths  are  given  by 
o l 


ray:  (z 


-2*1 


cos  0, 


sin  x 


Qq  ray  : same  as  equation  (8) 
where 
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V cos  3^ 


-V 

o 

V'  cos  9 

c 


We  wish  to  compute  Az  z - z for-  identical  values  of  at.  Solutions  of  the 
equations  for  each  ray  are: 


From  these,  it  is  evident  that 


3z 

dO 


lim 


Az 

A9 


dz 

!o 

5 

1+2  tan  0 I 

V'x  \ 

r \a  ] 

dO 

0 

V' 

b 0 

0 

O 

K j 

l Vo  / J 

1/2 


wue re  6 denotes  differentiation  with  respect  to  explicit  variables  only 


dz 


dO 


= x esc  0 


V'x 


I V'x 


1+2  tan  0 


-1/2 


(9) 


Therefore,  equation  (l) , becomes 
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But . from  equation  (7b),  one  can  compute  that  the  radical  is  nothing  more 
than  cos  0 sec  0Q,  therefore, 

P cos  Q sec  9 

T C 

x ^ __ g 

x sec  t cos  0 
o 

p (10; 

I = 

2 2 n 
x sec  0 


Equations  (2)  and  (.10)  aie  identical  in  form.  However,  in  the  interpre- 
tation of  equation  (10),  the  distance  along  the  raypath,  R,  is  no  longer 
equal  to  x sec  0 . Consequently,  tbs  raypath  intensity  contours  differ 
from  the  IsovelcSity  case. 

What  has  beers  derived  to  this  point,  equation  (10),  is  net  new;  it  Is, 
for  example,  derived  more  elegantly  in  C.  B.  Officer’s  "Introduction  to 
the  Theory7  of  Sound  Transmission,"  McGraw-Hill,  1958,  pgs-  48-50,  59 -6l. 

Continuous  Linear  Segment  Approximation  to  the  Non-linear  Gradient • 

It  has  been  customary  for  sotae  time  to  approximate  non-linear  sound-  v 

velocity  gradients  by  continuous  linear  segments  for  the  purposes  of  approxi- 
mate calculation  of  sound  raypaths  in  the  ocean.  The  ray paths  which  resuLt 
are  arcs  of  circles,  as  above,  but  with  different  radii  in  each  depth  seg- 
ment. The  tangent  lines  to  the  resulting  raypath  are  continuous  over  the 
boundaries  of  the  depth  segments.  In  general,  the  linear  segment  approxi- 
mation is  an  "eye"  fit  with  unequal  depth  segments. 
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The  equations  governing  the  segments  ere: 


V(z±) 

II 

< 

H- 

7{z) 

« \ + (*  * 

O-i) 

7i+l 

= v.  , <z.  + x -z±)  V' 

vhere  ail  subscripts  ran  over  i = 0,  1, . . . and  where  the  second  of 
equation  (11)  applies  only  between  z..  and  z.  + . 

dz/ud  for  Continuous  Linear  Segment  Approx  irsati  on 
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Here,  again,  ye  investigate  the  Unit  of  A*/A8.  but  for  the  situation 
of  the  raypath  generated  by  the  linear  uagaent  apprSxLnat lor. . The  notation 
i3  different  frca  the  above;  the  mania*?  subscript  n refers  to  the  n-ih 
depth  ssgaent  crossed  by  the  ray  vtalle  zr  re  ’era  to  the  initial  depth  of 
the  ray  (transducer  depth)  which  is  not  hccea3erily  or  a segment  bounuar^ . 
The  ray  picture  aav  look  c.ccte thing  like  Figure  3* 


In  Figure  3>  the  boundaries  of  the  depth  segment,  are  shown  by 
dashed  lines  intersecting  the  z -axis . We  observe  that  the  sequence  x 
is  increasing  while  the  sequences  z^  and  ©a  are  not  necessarily  so.  E 
If,  as  is  drawn  in  Figuie  3,  the  raypath  loaches  a point  of  tangeney  with 
a oegaant  boundary,  the  ray  la  considered  not  to  have  crossed  the  boundary 
and  will  continue  on  in  the  same  velocity  gradient  region  aa  before  the 
point  of  tangency. 

For  the  calculation  of  the  ratio  Az / AG  , the  notation  of  Figure  4 
applies . 0 


Geometry  and  Rotation  for  Az/AQ  Calculation 

o 


Figure 
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Care  nrust  be  taken  in  using  the  geometry  of  Figure  k because  of  the 

different  relative  magnitudes  of  x and  x! . Tie  rule  I3 

n n 

0 (positive)-*  x'  - x < 0 (ns  showu) 
n n n 

(12} 

9q  (negative  )->-  - xq  > 0 (not  ohovn) 

From  the  figure,  the  equations  governing  the  arcs  of  the  9 -ray (and  9 ray} 
frean  xq  (and  x^)  to  x are  given  by  equation  (13),  (respectively) 

2 V 2 V 

{z  - z )*  +(x  - X )2  + — (2  _ z ) E tan  0 (x  - x ) = 0 

Q y t C ^ ^ 

u a 


2 y 2 V 

(z  - zn)2  + (x  - X^)2  + (2  . z^) E tan  6n  (x  - x^)  * 0 

Vn  Vn 


From  Snell's  lav,  along  each  curve  one  gets 


n 


where  the  sign  is  positive  for  0 positive  end  negative  for  © negative.  The 

n n 

solution  of  equation  (13)  for  2 and  2 yield : 
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By  noting  that  x - x ' = >:  - x + Ax  with  the  same  convention  as  is  in 
(12)  and  then  expand?ng  the  fSrat  of  the  radicals  as  a function  of  Ax, 
one  gets 


Az 


A9  A9 

o o 


*<V  - X"V  j ♦ 


x(50> 


fa 


2 x „2 
sec  9 - V 
o n 


- — (x  - X ) 


Ax 


+ (Ax)  (ate) 


where  the  definition  of  X(9q)  is 


1 . 


2v;  fa 


sec2  0 - V2 

o n 


(X  - xn) 


fv;  (x 


xn^ 


1/2 


If  we  multiply  through  by  one  over  A0q  and  pass  to  the  limit,  we  arrive  at 
the  preliminary  form  of  the  recursion  1 onmt-a 


dx  Vn  i X(0o) 


dGo  Vn 

o n 


*9 


X(60) 


fa 


sec2  0 - V2 

o on 


_ V'  (x  - x) 

; -? L. 


dx 


d©_ 
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The  procedure  is  now  straight  f oruard : compute  the  partial  derivative, 
express  the  derivative  dz/dG(  in  terms  of  dsQ/d©Q,  and  simplify 


5X 

Wo 


V*  V2  sec2  G tan  9 (x  - x ) 
no  o o a 

Vn  X(G  ) V \/v2  sec2  9 - V2 

o n y o 0 n 


The  expression  for  dz  /do  can  be  evaluated  as  follows : 

no 


(15) 


Az  2 - z 

= for  0 positive  and 

A©  A©.  R 

o o 


Az  z - z 

— Fi  _ L_  for  © negative,  ’.there  the  aquations  for  z and  z are: 

A©  A©_  n 

o c 


2 v 2 V 

( z - z )2  + (Ax)2  + — (z  - z ) + — ten  © (Ax)  = 0 

33  Bl  a Vt  11 

n - 1 vn  - 1 


2 V 2 V 

(z  - z )2  + (Ax)""  + R — ( z - z„)  + — tan  ©^  (Ax)  = 0 


V* 

vn  - 1 


33  V'  n 

n - 1 


It  follows  then  that  in  the  limit  as  A©c  goes  to  zero,  the  formal 
expressions  for  the  derivatives  will  be  the  negative  of  each  other.  There- 
fore, the  Joint  calculation  can  be  written  as 


tan  ©n  (Ax) 
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Where  the  sign  convention  is  the  upper  when  is  positive  and  the  lower 
when  ©a  is  negative.  Expanding  the  radical  ai  a function  of  Ax  yields 


_ v»  j 

-"--1  tan  9 l Ax  + (Ax)2  (etc.) 

‘I 


1 - 1 


which  upon  passing  to  the  limit  becomes , exactly. 


= + ten  © 


or  = I tan  9 


in  which  the  sign  convention  is  the  usual  one.  This  is  certainly  the 
expected  result  for  this  calculation.  Using  equations  (l^),  (15),  (16),  and 


the  definition  of  tan  0^,  one  gets 


d#c  X(B0) 


i . VA  U ~ *t)  % 

V_  tan  9 d© 
a n J o 


tan  9 (x  - x ) 
o n 

cos  9 sin  9_ 
n n 


This  is  the  recursion  relationship  which  carries  the  derivative  over  the 
boundary  from  one  gradient  to  another.  To  further  simplify  the  equation, 
the  following  two  relationships  are  noted: 

X(0Q)  = cos  9 see  ©o 


Vn  *x  " xn*  #ln  9 


V„  tan  ©n 
n n 


sin  © 
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and,  therefore. 


or 


dz  cos  9n 

j sin  6 dz 

1 n 

tan 

d©  cos  0 

O 

sin  © d© 
n 0 

k 

cos 

dz 

dz 

sin  9 

M/>Q  f1  - _ fm 

oa  q El 

COb  11  1 ’ s C 

dG 

UD 

° d9 

sin  9 

0 

0 

L n 

^ eo  <x  - xn) 
sin  9_ 


(18) 


Repeated  application  of  equation  (l8),  decreasing  the  index  n,  yields  the 
desired  f crania : 


dz 


cos  0 


= tern  © 


dO. 


I 


(xj  - xi  _ L)  sin  0 


■f  sin  ^ sin  0^^ 

The  reciprocal  of  the  intensity  function  is  given  by 


x - x_ 


sin  0 


(19) 


1 x sin  0 

o 

I P cos  & 

o 

Equation  (20)  lias  the  following  properties; 

I.  It  reduces  to  the  inverse  of  equation  (10)  for  the  first  tern  only  case: 

, 2 
1 x 

— = g 
I P cos  9 
o 

II.  The  Bua  of  tvo  adjacent  terms  of  the  series  are  identically  equal  to  a 
corresponding  single  term  if  and  only  1£  the  velocity  gradients  across 
the  boundary  are  the  same : 


V"*  ^Xi  ‘ xi  - 1*  3ln  6 
Y*  sin  sin  9 ^ 


x - 


sin  9 


(20) 
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xn  + 1 ~ xn  , xn  - xn  - 1 . xn  + 1 ~ xn  - 1 

sin  9 . , sin  & sin  0 sin  0 , sin  9 , ain  9 , 

n+1  n n n-1  n+1  n-1 


if  and  only  if  V'  ==  V'  . 

n n-l 


III.  The  inverse  intensity  can  be  written  as  a sum  of  a tern  which  reduces 
to  the  linear  gradient  loss  plus  a sum  of  gradient  dependent  terms 
which  becomes  zero  in  the  absence  of  a change  of  gradient : 


1 

I 


2 

P cos  9 


- VQ  (sin  9 - sin  «0  ) 

V’  cos  9 
o o 


+ Ein  0,  sin  9 
o 


I 


sin  Qj  sin  9 

X 


(21) 


where  x_  . - x and  9 , = 0.  Thau  the  first  term  reduces  to  x when 

n+i  n+l 

the  gradient  Is  linear  can  be  seen  frcm  (7b). 

Equation  (20)  is  not  really  new  either;  its  essence  can  be  found  in 
equation  (89)  of  Chapter  3 of  "Physics  of  Sound  in  the  Sea"  and  elsewhere. 
However,  I have  not  seen  it  written  in  this  form  --  usually  it  is  written 
from  crossing  to  crossing  of  the  velocity  break  points  --  nor  have  I seen 
the  condition  that  only  actual  crossings  of  break  points  (and  not  tangencies 
to  them)  need  be  considered  in  the  summation . 
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